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ABSTRACT 

We test the concordance ACDM cosmology by comparing predictions for the mean 
properties of galaxy clusters to observations. We use high-resolution iV-body sim- 
ulations of cosmic structure formation and semi-analytic models (SAMs) of galaxy 
formation to compute the abundance, mean density profile, and mass of galaxy clus- 
ters as a function of richness and luminosity, and we compare these predictions to 
observations of clusters in the Sloan Digital Sky Survey (SDSS) maxBCG catalogue. 
We discuss the scatter in the mass-richness relation, the reconstruction of the cluster 
mass function from the mass-richness relation, and fits to the weak-lensing cluster 
mass profiles. The impact of cosmological parameters on the predictions is investi- 
gated by comparing results from galaxy models based on the Millennium Simulation 
(MS) and another WMAP1 simulation to those from a WMAP3 simulation. We find 
that the simulated weak-lensing mass profiles and the observed profiles of the SDSS 
maxBCG clusters agree well in shape and amplitude. The mass-richness relations in 
the simulations are close to the observed relation, with differences < 30%. The MS 
and WMAP1 simulations yield cluster abundances similar to those observed, whereas 
abundances in the WMAP3 simulation are 2-3 times lower. The differences in cluster 
abundance, mass, and density amplitude between the simulations and the observations 
can be attributed to differences in the underlying cosmological parameters, in partic- 
ular the power spectrum normalisation ag . Better agreement between predictions and 
observations should be reached with a normalisation 0.722 < erg < 0.9 (probably closer 
to the upper value), i.e. between the values underlying the two simulation sets. 

Key words: galaxies: general - galaxies: clusters: general - large-scale structure of 
the Universe - gravitational lensing - cosmology: theory - methods: numerical 



1 INTRODUCTION 

Clusters of galaxies are a powerful probe of astrophysics and 
cosmology (Voit 2005). For example, the cluster mass func- 
tion is very sensitive to the cosmic mean matter density, to 



the initial fluctuation amplitude ( 


Press and Schechter|1974| 


Frenk et al.|1990 Eke et al.||1996 


Sheth and Tormen|1999|) 


and to dark energy dynamics ( Bartelmann et al. 2006 Grossi 


and Springel|2009| Francis et al.||2009|). It can be predicted 



to high accuracy by numerical simulations ( jJenkins et al 
2001||Warren et al.|2006| |Lukic et al.|2007[ ). 

However, the masses of galaxy clusters cannot be ob- 
served directly. Thus one either needs to infer the masses of 
observed clusters from some more directly observable cluster 



property, requiring accurate knowledge of the observable- 
mass relation and its scatter, or one must directly compare 
observed and predicted cluster abundance as a function of 
such observables. These include the X-ray luminosity and 
temperature of the intracluster gas (Borgani et al.||2001| 



Reiprich and Bohring 


er 2002; Stanck et al. 2006, Piffaretti 


and Valdarnini 


2008 


Vikhlinin et al. 


2009|, the number 



and velocity dispersion of the cluster ga laxies (|Zwicky|1937| 



Rines et al.|2003||KocIianek et al.|2003||Becker et al.|2007j, 



the number of giant arcs ( Bartelmann et al.|1998 Wambs- 
ganss et al. 2004; Fedcli et al. 2008), the weak-lensing signal 



(Tyson et al. 1990 Cypriano et al. 2004 Hoekstra 2007 



Johnston et al.|2007b Reyes et al.|2008 1 , and the Sunyaev- 



shilbertSastro . uni-bonn . de 



Zel'dovich signal induced by the cluster ( White et al.||2002 



Schulz and White||2003l |Bonaldi et aLl|2007| |Staniszewski| 
et al.|2009| |~ 
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Predictions for many cluster observables (e.g. the X- 
ray luminosity or the cluster richness) and for their relation 
to the cluster mass require modelling of astrophysical pro- 
cesses such as gas cooling and galaxy formation. Although 
inaccuracies in such astrophysical models are an unpleasant 
source of uncertainty for cosmological parameter estimation, 
they mean that comprehensive observations of clusters can 
be used to constrain cosmological parameters and models 
for cluster/galaxy evolution simultaneously. 

The largest sample of observed galaxy clusters currently 



available is the maxBCG cluster catalogue (Koester et al 



2007a), This was extracted from the Sloan Digital Sky Sur- 
vey (SDS S)p| using max BCG, an optical cluster-finding al- 



gorithm (Koester et al. 2007b I. Constraints on the scatter 



in the velocity dispersion-richness relation and the mass- 
richness relation of these clusters have been derived from 
cluster X-ray and galaxy velocity dispersion observations 
( Becker et aLl|2007| [Rozo et al.||2009a[ ). Weak-lensing mea- 
surements of average cluster mass profiles have also been 
used to calibrate the mass-richness relation, the mass-optical 
luminosity relation, and mass-to-light ratio profiles of the 
maxB CG clusters (|Sheldon et al.||2009a|b| | Johnston et al.| 



2007b Reyes et al. 20081. These data provide significant 



constraints on cosmological parameters ( Rozo et al.|20 09b) 



In this work, we investigate how well physically based 
models for galaxy formation in a ACDM universe can re- 
produce the observed relations of cluster richness and lumi- 
nosity to other cluster properties, most notably mass. We 
also investigate what information on cosmological param- 
eters and galaxy evolution can be obtained by comparing 
model predictions to observation. We use the Millennium 



Simulation ( Springel et al. 2005 1 and two smaller iV-body 



simulations of cosmic structure formation ( Wang et al. 2008 1 
in conjunction with semi-analytic models of galaxy evolu- 
tion ( |De Lucia and Blaizotl|2007| |Wang et al~1|2008[ ) to cre- 
ate mock catalogues of galaxy clusters selected similarly to 
the maxBCG catalogue. We compute cluster abundances, 
average cluster masses, and weak-lensing mass profiles as a 
function of cluster richness and luminosity, and we compare 
these to observational results for the SDSS maxBCG sample. 
In addition, we investigate the scatter in the mass-richness 
relation, and we discuss how well one can recover the cluster 
mass function and the weak-lensing mass profiles from the 
richness-binned cluster abundances and mean masses. 

The semi-analytic galaxy models used here couple star 
formation in galaxies to the properties of the evolving dark 
matter halo distribution in which the galaxies live. The mod- 
els have been adjusted to be consistent with various obser- 
vations, e.g. the luminosities, stellar masses, morphologies, 
gas contents and correlations of galaxies at low redshift, but 
they have not been tuned to match the properties of rich 
clusters. The comparison to observations provided here is 
thus a direct test of these models and their description of 
the physical processes relevant for galaxy formation. This 



contrasts with halo occupation distribution models ( Cooray 
|and S hcth 2002]) , where the galaxy populations of clusters 
are adjusted to fit observation without considering in detail 
how they could have been built up by physical processes 
within the evolving dark matter distribution. 



Table 1 . The cosmological parameters at redshift for the three 
simulations used in this study. The parameters are: the baryon 
density £7^, the matter density SImi and the energy density of the 
cosmological constant Qjy (in units of the critical density), the 
Hubble constant h (in units of lOOkrns - ^^Mpc -1 ), the primordial 
spectral index n and the normalisation parameter ag for the linear 
density power spectrum. 





MS & WMAP1 


WMAP3 




0.045 


0.04 




0.25 


0.226 


n A 


0.75 


0.774 


h 


0.73 


0.743 


n 


1 


0.947 


OS. 


0.9 


0.722 



Our paper is organised as follows. We discuss the N- 
body simulations and galaxy models, as well as our meth- 
ods for creating the simulated cluster samples from them in 
Sec. [2] Results for our simulated cluster samples and a com- 
parison to observation are presented in Sec. [3] Our paper 
concludes with a summary and discussion in Sec. [4] 



2 METHODS 

We use cosmological iV-body simulations to analyse the mat- 
ter distribution in and around galaxy clusters in two differ- 
ent ACDM cosmologies. We infer the properties of the galax- 
ies in the clusters from model galaxy catalogues generated 
by applying semi-analytic galaxy formation models to halo 
assembly trees generated from the outputs of the A r -body 
simulations. We then compute richness and luminosity es- 
timates for clusters in the model galaxy catalogues taking 
into account several observational features of optical cluster- 
finding algorithms, in particular maxBCG by Koester et al. 
(2007bt. 



2.1 The A-body simulations 

Our study is based on three different A r -body simulations: 
the Millennium Simulation (MS) by |Springel et~aL] ( |2005[1, 
and two sma ller simulations WMAP1 and WMAP3 by |Wang| 



et al 



_(2OO8)0The simulations assume a flat ACDM cosmol- 
ogy with parameters given in Table [I] Both the MS and the 
WMAP1 simulation use a parameter set that was derived by 



combining the WMAP lst-year results ( Spcrgcl et al.|2003 ) 



with results from the 2dFGRS ( |Colless et al.| |2003). The 
WMAP3 simulation employs cosmological parameters that 
are consistent with data from the WMAP 3rd-year release, 
the Cosmic Background Imager, and the Very Small Array 
(Spergel et al. 20071, with a bias towards values differing 



from those used for the other two simulations. 

The most prominent differences between the two sets of 
cosmological parameters are the normalisation parameter as 
and the spectral index n for the density power spectrum: The 
MS and the WMAP1 simulation assume as — 0.9 and n = 1, 
whereas the WMAP3 simulation assumes lower values as = 



1 http://www.sdss.org 



2 We refer the reader to |Springel et"aL] | |2005| , |Wang et al.1|2008| l, 

and references therein for more details about the simulations. 
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Table 2. Numerical parameters for the three simulations used 
in this study. The parameters are: the comoving cube size L, the 
particle number n p , the particle mass m p , and the effective force 
softening length e. 





MS 


WMAP1 


WMAP3 


L [/i^Mpc] 


500 


125 


125 


n p 


2160 3 


560 3 


560 3 




8.6 x 10 8 


8.6 x 10 s 


7.8 x 10 s 


e [h-^pc] 


5 


5 


5 



0.722 and n — 0.94. Hence, there is less power on small scales 
in the matter power spectrum of the WMAP3 simulation 
than in the MS and the WMAP1 simulation. This results in 
a substantial delay in structure formation and less massive 
collapsed structures at any given redshift in the WMAP3 
simulation. 

The simulations were run using a parallel TreePM ver- 
sion of GADGET2 ( |Springel|20"05t . The numerical param- 
eters of the simulations are listed in Tabled The main dif- 
ference between the MS and the WMAP simulations is the 
simulation box size. The large volume of the MS provides us 
with a large sample of galaxies and galaxy clusters with sta- 
tistical errors comparable to the 'cosmic variance' errors of 
the SDSS maxBCG sample. The smaller WMAP simulations 
differ in their cosmological parameters, but share the same 
numerical parameters and initial conditions]^] This reduces 
the influence of sampling noise when comparing results be- 
tween the WMAP simulations and allows us to study the 
influence of cosmology on our results. 

For each simulation, the particle data were stored on 
disk at 64 output times. These snapshots contain informa- 
tion on dark-matter halos, which have been identified run- 
ning a friend-of-friend (FOF) group finding algorithm on 
the set of simulation particles. The halos were later decom- 
posed into subhalos using SUBFIND ( |Springel et al.||2001 1 
to identify gravitationally self-bound locally overdense re- 
gions. The most massive subhalo, called the main subhalo, 
typically contains 90% of the FOF halo mass and shares its 
centre. Detailed merging history trees of all self-bound dark- 
matter subhalos were then computed. The resulting merger 
trees were used as input for the semi-analytic models dis- 
cussed in the following section. 



2.2 The semi-analytic galaxy models 

We use semi-analytic galaxy formation models from the Mu- 



nich family 


Kauffmann et al.|1999 Springel et al.|2001 


De 


Lucia et al.| 


2004| Springel et al.||2005| |Croton et al.||2006| 


De Lucia and Blaizot| 20071 to set the optical properties 



of galaxies in the A?-body simulations. These semi-analytic 
models assume that galaxies form from gas that accumulates 
in the centre of each dark-matter subhalo in the simulation. 
Star formation in each galaxy is coupled to its subhalo prop- 
erties via simple prescriptions for gas cooling, star formation, 



3 Their initial density fields are identical except for small ampli- 
tude adjustments needed to reproduce the correct matter power 
spectra. 



chemical enrichment and feedback from supernova and cen- 
tral black holes (AGNs), as well as for merging of galaxies 
once their dark matter subhalos have merged. Certain pa- 
rameters quantifying the efficiency of these processes can be 
adjusted in order to maximise agreement of the results with 
observation. 

For the MS, we use the publicly available galaxy cata- 
lo gu^]that was generate d using the galaxy model described 
in De Lucia and Blaizotj ( |2007[ ). This MS model reproduces 
various observed relations for galaxies, in particular the ob- 
served luminosity, colour, gas content and morphology dis- 



tributions (Croton et al. 2006 De Lucia et al. 2006 Kitzbich 



|ler and White|2007|l and the o bserved two-point correlation 
functions ( |Springel et al.|2005||Kitzbichler and White|2"007 |. 

The galaxy model of De Lucia and Blaizot (20071 has 
been applied by |Wang et al.| ( |2008[ > to the WMAP1 simu- 
lation using the same set of efficiency parameters. As ex- 
pected, this gives similarly good agreement with observa- 
tions as for the MS. Here, we will use these model galaxies, 
which we refer to below as the WMAP 1- A model, for the 
W MAP1 simulat i on. 

|Wang et al.| (|2008[> also applie d the galaxy modelling 
technique of |De Lucia and Blaizotl (120071) to the WMAP3 



simulation, but with two slightly different sets of efficiency 
parameters. One, which we call WMAP3-B here, employs 
the same star formation efficiency as the WMAP1-A model, 
but lower supernova and AGN feedback efficiencies. The 
other model, called WMAP3-C in the following, employs 
a higher star-formation efficiency but also higher feedback 
efficiencies. Both the WMAP3-B and the WMAP3-C model 
show good agreement with low-redshift galaxy observations, 
e.g. of the galaxy luminosity function and the galaxy clus- 
tering ( Wang et al.|2 008). In the following, we will consider 
both models for the WMAP3 simulation. 



2.3 The ridgeline galaxies 

Optical cluster finding algorithms such as maxBCG ( jKoester] 



et al. 2007b I locate galaxy clusters by searching for local 
overdensities of E/S0 ridgeline galaxies in angular and red- 
shift space. These galaxies are common in known clusters, 
and they are relatively easy to find, since many of them 
are bright and their observed colours are strongly correlated 
with luminosity and redshift. Because of their small scatter 
in colour, a narrow search range in colour can be chosen at 
each redshift. This narrow search range forces one to know 
the mean ridgeline colours rather accurately. Accurate mean 
ridgeline colours are also important for accurately quantify- 
ing the ridgeline galaxy content of the clusters. 

Although the semi-analytic models can reproduce many 
of the observed properties of galaxies, there are still some 
discrepancies. In particular, the models do not reproduce 
the colours of passively evolving galaxies to the degree re- 
quired for a direct application of the ridgeline colour-redshift 
relation used for the maxBCG catalogue (see, e.g 



Wein- 



|mann et al.||2006] for possible reasons). We therefore 'mea- 
sure' the mean ridgeline colours of model galaxies as a func- 
tion of redshift: We roughly identify the mean colour of the 
ridgeline galaxies in the colour-magnitude diagram (where 



http://www.mpa-garching.mpg.de/Millennium 
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ilaxies as a 



Figure 1. The mean colour g — r(z) of the ridgeline j 
function of redshift z for the galaxy model of |De Lucia and Blaizot| 
l |2007| in the MS compared to the mean colour-redshift relation 
for the ridgeline galaxies in the SDSS (Ko ester et al.| 2007b). 



the ridgeline population induces a visible overdensity among 
the bright galaxies with colours close to the observed mean 
ridgeline colour). We then fit a Gaussian with mean x and 
variance a 2 to the distribution of galaxy colours in a region 
around our initial colour guess. We repeat the fit considering 
all galaxies with absolute i-band magnitudes Mi ^ —20 and 
colours in the range x ± 3a until x converges]^] 

The resulting mean ridgeline colours for model galaxies 
in the MS as a function of redshift is compared to the cor- 
responding relation for SDSS galaxies in Fig. [I] The mean 
ridgeline colours of the model are close the observed ridge- 
line colours at low redshift, but they deviate significantly at 



higher redshift. The measured ridgeline width a g 



0.05 



for all considered galaxy models and redshifts, which is in 



agreement with observations (Koester et al. 2007b I. The 
measured width in (r — i)-colour ov_; « 0.03, which is 
smaller than the observed value a r -i = 0.06. 

Following the maxBCG observational procedure, we 
consider g — r and r — i colours to identify ridgeline galaxies 
in our simulations. For each simulation snapshot with red- 
shift 0.1 ^ z ^ 0.3, we select all objects in our semi-analytic 
galaxy catalogues whose g — r and r — i values are both 
within 2a of the mean ridgeline colours. As mean colours, 
we take the values measured from the simulations. For the 
ridgeline width a, we take the observed values a g - r = 0.05 
and a r -i = 0.06 ( |Koester et aL|2007b |. 

Besides colour, galaxy brightness is used to select ridge- 
line galaxies. We thus further select from all model galaxies 
surviving the colour selection those with apparent observer 



„obs 
*lim 



$ j Here, we employ the 



Koester et al. 



(2007b I (B 



frame i-band magnitude i 
same magnitude limit 
Koester, private communication). This magnitude limit cor- 
responds to an absolute rest-frame magnitude limit A/, rcst « 
—20.25 + 5 log 10 h « —20.9 for the cosmologies considered 
here. 



2.4 The galaxy clusters 

We identify galaxy clusters in the simulations with dark mat- 
ter halos found by the FOF algorithm. From the simulation 
data stored on disk, we obtain for each such halo (hence 
cluster candidate) the positions of the centre, the virial ra- 



dius R 



.crit /mean 
200 



(i.e. the radius of the largest sphere enclosing 



200 x the critical/mean density) and the mass Mggo ' mean 
(i.e. the mass enclosed within R^^ meaXL ). 

The semi-analytic galaxy models provide us with in- 
formation about the galaxies associated with each dark 
matter halo. For each halo, we measure the total number 
Nfa t of associated ridgeline galaxies selected as described 
in the preceding section. In addition, we count the num- 
ber of ridgeline galaxies Nf^ within a physical radius of 
lh -1 Mpc in projection along a simulation box axis. The re- 
sult is used to compute the "observatio nally defined" radius 



flgft = 0T56(iySL) - 6 fc~ X Mpc (s 



Hansen et al. 



2005 



Koester et al. 2007b I. We then calculate a scaled galaxy 



richness N^q Q by counting all ridgeline galaxies within a pro- 
jected radius -Rfoo- A cluster luminosity Lfao is then com- 
puted by summing the i-band luminosities of all ridgeline 
galaxies within Jifoo- These procedures mimic closely those 
employed to estimate richnesses and luminosities for the real 
maxBCG catalogue. 

To improve statistics, we perform the measurements of 
the projected quantities Nf M , R^qq, -^foo> and L|q along 
all three simulation box axes. Each projection is treated as 
an individual cluster in the subsequent analysis. For the 
computation of cluster densities, the resulting triplication 
of clusters is taken into account by assuming a three times 
larger simulation volume. 

Several effects hamper a direct comparison between ob- 
served clusters and our simulated clusters at the stage de- 
scribed so far. Here, we will not take into account fragmen- 
tation, but we do correct for contamination of the observa- 
tional data by foreground and background galaxies)^] Using 
spectroscopic data, Koester et al. ( 2007a I found that roughly 



16% of the galaxies identified by maxBCG as cluster ridge- 
line galaxies are, in fact, projections. We include such pro- 
jections in our simulated clusters in a very simple way, by 
randomly duplicating about 19% of the ridgeline galaxies. 
As a result, our model clusters appear to be contaminated 
at the 16% level. 

We also consider the effect of misidentifying the clus- 
ter centre. For each cluster in the simulations, we calculate 
galaxy numbers, radii, etc. not only using the 'true' cen- 
tre, but also using the position of the second-most massive 
subhalo as 'apparent' centre. The resulting distribution of 
projected offsets R 3 between 'true' and 'apparent' centre is 
shown in Fig.[2]for the MS. The distribution can be approx- 
imated by a two-dimensional Gaussian distribution with 



pdf(-R s ) = ^ exp 



2a 2 



(1) 



and a s = 0.38h Mpc, which agrees well with the offset 



5 The dependence of the mean ridgeline colour on magnitude is 
ignored in this procedure, since for all models the slope in the 
colour-magnitude relation is very small in the considered colours 
and magnitude range. 



6 Fragmentation is insignificant in the maxBCG sample, but 
overmerging slightly boosts the cluster richness estimates mainly 
due to contamination by foreground and background structures 
I Koester et al.|2007b l. 
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Figure 2. The distribution of projected offset R s between the 
'true' and 'apparent' centres of miscentred clusters in the MS 
(dashed line). The distribution is well fit by a 2D Gaussian (solid 
line). 



distribution found for the maxBCG algorithm run on simu- 



lated data ( Johnston et al.||2007b |. We find a similar value, 
a s — 0.41/i~ Mpc, for the WMAP1 simulation, and some- 
what smaller values, a s = 0.34h _1 Mpc, for the WMAP3 
simulation. 

When needed for our analysis, we will use a probability 



1 + 0.04JV/ 
2.2 + 0.05N! 



gal 



gal 



(2) 



that a cluster in the simulation with ridgeline galax- 

ies is correctly centred. Empirically, this yields roughly the 
same probability PciN^o) that a cluster with measured rich- 
ness A^|qq is correctly centred as was found for the maxBCG 
algorithm by |Johnston et al. ( 2007b I . 

In our simulations, centre misidentification tends to re- 
duce the number iVfoo of ridgeline galaxies within i?|oo f° r 
a given cluster. Consequently, the number density of clus- 
ters with JVIqq above a given threshold is slightly decreased 
by centre misidentification. Another consequence of centre 
misidentification in our simulations is a slightly higher av- 
erage cluster mass and ridgeline galaxy number for 
a given measured richness N^qq- All these effects may be 
smaller for the actual maxBCG algorithm. This algorithm 
disfavours identifying the cluster centre with galaxies that 
lead to low ATfoo in comparison to galaxies that yield a larger 
cluster richness. In the following, we will thus discuss results 
for our simulated cluster samples in the case that centre 
misidentification is ignored, unless stated otherwise. 



2.5 The cluster samples 

A comparison of cluster abundances in our simulated cluster 
samples to observation requires knowledge of the volumes 
and areas of the real surveys. For the SDSS maxBCG cluster 
sample, we assume an effective survey area of 7400 de g 2 and 
a redshift range ofO.lsCz^O.3 ( |Rozo et al."||2009b[ ). This 
yields an effective survey volume of 4.3 x 10 a h~' i Mpc 3 for 
the MS and WMAP1 cosmology, and 4.4 x 10 s h~ 3 Mpc 3 for 
the WMAP3 cosmology. 

For each snapshot of our simulations with redshift 0.1 ^ 
z ^ 0.3, we create a model cluster catalogue containing the 



projections of all clusters in the simulation box. We calculate 
the statistical properties of interest (e.g. the cluster densities 
or average cluster mass as a function of cluster richness) 
for each of these snapshot catalogues (properly taking into 
account the increased cluster number due to inclusion of 
multiple projections of each cluster). We then compute the 
properties of a cluster sample in a volume-limited survey 
with 0.1 ^ 2 ^ 0.3 by an average over the snapshot results, 
where each snapshot is weighted by its cluster abundance 
and its volume fraction in the survey. 

The simulation cube volume L 3 = 1.25 x 10 8 h~ 3 Mpc 3 
for the MS, which is about one quarter of the SDSS maxBCG 
survey volume. The simulation cube is used in three differ- 
ent projections and at six different redshifts to construct the 
model cluster samples, which increases the effective sample 
sizeQOne can thus expect the statistical errors due to sam- 
ple variance to be roughly the same for the MS model and 
the SDSS maxBCG cluster sample, and about eight times 
larger for the WMAP models (with their 64 times smaller 
box volume) . 

Where appropriate, we estimate the errors of our models 
due to sample variance in the following way: We divide the 
simulation cube of the MS into 64 smaller cubes each having 
the size of the WMAP simulations. We calculate the observ- 
ables for each of these subcubes separately. The standard de- 
viation of the results from the different subcubes serves as 
an estimate for the statistical error of the WMAP 1- A model. 
The statistical errors for the other models are than extrap- 
olated from the WMAP 1- A error using simple assumptions 
about the scaling with volume, numbers, etc. 



3 RESULTS 

Here, we compare the properties of clusters in our various 
galaxy formation models to the observed properties of clus- 
ters in the SDSS maxBCG catalogue. We investigate average 
cluster properties as a function of galaxy content by dividing 
our simulated cluster samples into bins of richness iVfoo and 
luminosity Lfoo ( as was done for the maxBCG cluster sam- 
ple) . The average properties of the model clusters binned by 
Affgo are listed for the different galaxy formation models in 
Tables |7|10| The properties of mod el clusters when binned 



by Lfoo are shown in Tables ll]l4 



We first present results for stacked weak-lensing mass 
profiles and for the abundance of clusters, both as a function 
of galaxy richness. We then discuss the mean and scatter of 
cluster mass as functions of richness and luminosity. Finally, 
we study how well the cluster mass function and the stacked 
weak-lensing mass profiles can reconstructed from the clus- 
ter abundance together with the mean and scatter of the 
mass-richness relation. 



3.1 Cluster mass profiles 

For each cluster in the simulations, we compute weak-lensing 
mass profiles by projecting the simulation particles and the 



7 The resulting effective sample size is difficult to quantify since 
the subsamples created from the different projections and red- 



shifts are not independent. 
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galaxies in a cuboid region of 35/i _1 Mpc transverse physical 
side length and 100/i _1 Mpc comoving thickness centred on 
the cluster. The projected particles and galaxies (assumed 
to contribute their stellar mass as points J] are binned in 
annuli to compute the circularly averaged surface mass den- 
sity E(i£) at radius R from the projected cluster centre, the 
mean enclosed surface mass density E(-R) inside R, and the 
weak-lensing mass profile 



AE(7?) = £(i?) - £(#). 



(3) 



The weak-lensing mass profile AE(i?) is proportional to 
the average tangential shear (7 t ) (R) around the projected 
cluster centre, and can therefore be measured using weak- 



lensing observations (Schneider et al. 20061. To increase 



signal-to-noise, the measured tangential shear may be av- 
eraged over a sample of clusters. The resulting shear signal 
can then be converted to an average mass profile for the 
observed cluster sample. 

In this section, the weak-lensing mass profiles of the 
simulated clusters in each snapshot are averaged in bins of 
Affoo' The average profiles from each snaphot are then com- 
bined with appropriate weights to obtain a mean profile for 
each iVfQQ-bin in a volume-limited survey. In doing this, we 
also take into account the effect of cluster centre misiden- 
tification, which has a considerable im pact on the average 
profiles at radii R < Ih^ 1 Mpc (see Sec. 3.7 1. 

In Fig. |3J the weak-lensing mass profiles for the simu- 
lated clusters in the MS and the WMAP models are com- 
pared to the measured profiles of maxBCG clusters in the 



SDSS ( |Sheldon et al.||2009a[ ). The simulated and observed 
profiles agree remarkably well in detailed shape and ampli- 
tude. This is strong evidence that our models provide a re- 
alistic description not only of the density profile and galaxy 
content of galaxy clusters, but also of the maxBCG cluster 
selection and richness measurement. 

Differences between the galaxy models and the obser- 
vations are small but noticeable. The simulated density pro- 
files of the MS tend to be above the observed profiles in the 
poorer clusters but are an excellent fit in the rich systems. 
The WMAP3 model profiles fit better in the poor clusters 
but are mostly below the observed profiles in richer clusters. 
For radii R < l/i _1 Mpc, the mass profiles of the WMAP3 
models are roughly 30% lower than in the MS/WMAP1 
models. This suggests that for given richness, the maxBCG 
clusters are on average slightly less massive than clusters 
in the MS/ WMAP 1 simulations but slightly more massive 
than those in the WMAP3 simulations. We will investigate 
this issue in more detail in Sec. 13.31 



3.2 Cluster abundance 



ch- 



Table[3] compares the area density of clusters with ri 
ness A^Iqq ^ 10 and redshift 0.1 ^ z ^ 0.3 in our simulated 



8 Although particles in the simulations represent the total mass 
in the simulated parts of the universe, we do not compensate for 
the additional stellar mass, since (i) the mass in stars is very 
small compared to the total mass in collapsed objects, and (ii) 
gas physics increases the dark-matter density in the inner part of 
the halos compared to collisionless simulations (e.g. |Barnes and| 
White|1984||Gnedin et al.|2004 i. 



Table 3. The area density ra >10 [in deg ] of clusters with rich- 
ness A^Iqq 10 and redshift 0.1 $C z $C 0.3 for the model clusters 
in the MS and WMAP simula tions and for the SDSS maxBCG 
clusters | |Koester et al.|2007a^ . Considered are the cases that all 
cluster centres are correctly identified (p c = 1), and that only a 
fraction of clusters given by Eq.[5]is correctly centred (p c < 1). 
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Figure 4. The comoving abundance nS,"") (z) of clusters with 

richness N^qq 10 as a function of redshift z in the SDSS 
maxBCG catalogue (Koester et al. 2007a) and in our simulated 
cluster catalogues based on the MS and the WMAP simulations. 



surveys to the observed abundance of maxBCG clusters in 
the SDSS. If we ignore misidentification of the cluster cen- 
tres, the MS and the WMAP1-A model yield cluster abun- 
dances very similar to the real maxBCG cluster sample. In 
contrast, the cluster abundances for the WMAP3 models 
are a factor 2-3 lower than that observed in the SDSS. The 
lowest abundance is found for the WMAP3-B model. 

To obtain an estimate for the statistical error on the 
cluster density, we employ the subsampling method de- 
scribed in Sec. 12.51 We use the MS model to create subsam- 
ples with sizes equal to the WMAP 1- A model. From these 
subsamples, we estimate a standard deviation for the area 
density of clusters of 25% for the WMAP1-A model. Tak- 
ing into account the lower cluster densities in the WMAP3 
models, a slightly larger statistical error of 30-40% can be 
assumed for these. A simple extrapolation of the WMAP1- 
A error yields a much smaller error of 3% for the cluster 
density in the MS model. 

Taking into account centre misidentification in our sim- 
ulations reduces the cluster abundances by « 20% and 
thus increases the discrepancy between the WMAP3 models 
and the SDSS cluster sample. The abundance is then also 
somewhat low even in the MS case, suggesting, as noted 
above, that our procedure for modelling the effects of cen- 
tre misidentification may overestimate the effect in the real 
maxBCG catalogues. 
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Figure 3. The average weak-lensing mass profile AE(iJ) as a function of radius R for clusters in different richness bins in the MS 
(dashed lines), WMAP1-A (dotted lines), WMAP3-B (dash-dotted lines), and the WMAP3-C model (dash-dot-dotted lines) compared 
to the observed profiles in the SDSS (She ldon et al.|2 009a points with error bars). 



8 



S. Hilbert & S.D.M. White 



CI j 

1) 



o 



10 

1 

0.1 



51 10" 



10 



-3 



10 



-4 



SDSS 

MS 

WMAP1 A 

WMAP3 B 

WMAP3 C 



10 1 



A' 



10 



20 



50 



, \ 

\ 
100 



N: 



gal 
2(ill 



Figure 5. The surface density n ang (A r |Q ( l ) ) of clusters with red- 



shift 0.1 $J z < 0.3 as a function of richness AT; 



Kill 
200' 



Counts in the 



SDSS maxBCG cluster sample jKoester et al.||2007a| ) are com- 
pared to counts in our simulated cluster catalogues for the MS 
and the WMAP simulations. 



The redshift dependence of the comoving cluster abun- 
dance is shown in Fig. [4] All our models show a slight de- 
crease of the comoving cluster density with increasing red- 
shift. The statistical errors on the comoving cluster abun- 
dance estimated with the subsampling method are similar 
to the errors on the cluster area density, i.e. 25-40% for 
the WMAP models and ~ 3% for the MS. For z > 0.15, 
the abundance in the MS and WMAP1-A model agree very 
well with the maxBCG results. In contrast, there is a much 
larger observed density at low redshifts z < 0.15, which is 
not seen in the models and presumably reflects nearby large- 
scale structure such as the SDSS Great Wall. 

The dependence of the cluster counts on richness Af 00 
is illustrated in Fig. [5] Cluster abundances in the MS and 
WMAP 1- A models are lower than in the SDSS for Ar| a < 



20, but exceed the observed abundances for A/foo ^ 20. 
Abundances in the WMAP3 models are always lower than 
in the MS and WMAP 1- A model and below the observa- 
tions. For Afoo > 9, the cluster abundances in the WMAP3 
models are 2-20 times lower than the abundances in the 
WMAP1-A model and 2-5 times lower than the observed 
abundances. Hence, the differences are always larger than 
the statistical errors inferred from the subsampling (which 
are ~10-100% for the WMAP models and ~1-10% for the 
MS and SDSS). The differences are largest in the high-A^fgQ 
tail of the distribution. 

The low cluster abundance for the WMAP3 models in 
comparison to the MS and WMAP 1- A models, in partic- 
ular for large A/fog, are a reflection of the different cos- 
mologies. Rich clusters have massive dark matter halos, and 
there are far fewer massive halos in the WMAP3 cosmology 
than in the WMAP1 cosmology. This is mainly due to the 
lower value of as- The higher star formation efficiency in the 
WMAP3 models does not sufficiently enhance the number 
of bright ridgeline galaxies per unit mass to compensate for 
the decrease in the number of massive halos. As a result, 



*s 10 1 



5 ri 

5 



10' 



















\ SDSS 






--MS 
--- - WMAP1 A 






- - - WMAP3 B 


£1 




■ - - WMAP3 C - 


5 


10 20 

<<> 


50 100 



Figure 6. The average cluster mass (A^foo) vs ' r i cnness (-^foo) 
relations for cluster catalogues from the MS and WMAP sim- 
ulations are compared to the relation derived from the SDSS 
maxBCG catalogues by |Johnston et al.| |2u07bj . 



there are fewer rich clusters in the WMAP3 models than in 
the MS and WMAP 1- A models. 

The observed abundance of rich clusters is slightly lower 
than the abundance predicted for the MS and much larger 
than the abundance predicted by the WMAP3 models. This 
suggests that 0.72 < as < 0.9 for our Universe, with as 
probably closer to 0.9 than to 0.72. This is consistent with 
some recent estimates [e.g. as 
[et~aLl[2007l <r 8 

as = (0.83 ± 0.03)(!T2m/0.25)- ' 41, 



: 0.80 ±0.02 by|Lesgourgues 
0.81 ± 0.03 by |Komatsu et al ||2009r and 
by |Rozo et al.|2009b 



The results for the abundances of rich clusters alone are 
not sufficient to definitely conclude that 0.72 < as < 0.9. 
For example, problems with the modelling of the ridgeline 
galaxies and their identification could have lead to inaccu- 
rate estimates for the cluster richness and thus to incorrect 
cluster abundances. However, there is complementary evi- 
dence from the mass-richness relation, which we discuss in 
Sec. rO and I3U 



3.3 Mass-richness relation 

The weak-lensing mass profiles discussed in Sec. |3.1| can be 
used to estimate spherical-overdensity cluster masses. This 
can be done, e.g., by a non-parametric conversion of the 
wcak-lcnsing mass profiles into average 3D density profiles 
( Johnston et al.|2007a l, or by fitting a parametrised model 
of the average cluster density to the shear data ( Johnston 
et al.|2007bl. 

In Fig. |6l average cluster masses (Affoo) m our simu- 
lated catalogue^] are shown as a function of cluster richness 
(A/Joo) and are compared to the corresponding relation for 



9 We use the masses measured directly from the matter dis- 
tribution in the simulations. As discussed in Sec. |3.7| these are 
consistent with the masses obtained from parametric fits to the 
weak-lensing mass profiles. 
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Figure 7. The average cluster mass 

( M 200 an > vs - richness (JVfgJ) 
relations for cluster catalogues from our MS and WMAP simu- 
lations are compared to the SDSS relation given by | Reyes et al,| 
( |2008| . 



SDSS maxBCG clusters by Johnston et al. ( 2007b I, who cal- 



culated the cluster masses by fitting parametric models to 
the observed lensing signalrj Remarkably, all our simula- 
tions reproduce the observed mass-richness relation within 
~ 30% over the two orders of magnitude spanned by the 
SDSS clusters. This is corroborates our finding in Sec. |3.1| 
that our models provide an adequate description of the sta- 
tistical properties of optically selected galaxy clusters. 

At given richness, the differences in cluster density pro- 
files between models and observations (see Fig. [3| imply 
differences in mean cluster mass. For ArfgJ < 10, the MS 
yields cluster masses that are 30-40% higher than the SDSS 



maxBCG cluster masses derived by Johnston et al. (2007b I. 
For ATfoo > 10, the MS masses are up to 20% higher than the 
SDSS maxBCG masses. The cluster masses of the WMAP1- 
A model are comparable to those of the MS, but seem to 
be affected by sampling noise for large (iVfoo). The cluster 
masses in the WMAP3-B model are lower than those for 
the MS model by 20-30%, and fall below the SDSS cluster 
masses at large richness. Among our models, the lowest clus- 
ter masses are found for the WMAP3-C model, where the 
values are always smaller than those in the SDSS. 

The mass-richness relation for the WMAP 1- A model 
differs by up to 30% from the relation for the MS model 
(which is much less affected by sampling noise due to its 64 
times larger volume) . Using subsamples created from the MS 
as described in Sec. 12.51 we estimate a standard deviation 



for the binned cluster mass of 2% for Nioo = 3 and 20% 
for 71 < Af a o < 220 for the WMAP 1- A model. A similar 
statistical error can be assumed for the WMAP3 models. 
This is consistent with the differences between the MS and 
WMAP 1- A models being solely due to sampling noise. 

Another analysis of the weak-lensing data for the SDSS 



10 We multiplied the masses given in | Johnston et ah] ( |2007b| | by 
a factor 1.18 ± 0.06, which is a photo- z bias correction advocated 



Table 4. The best-fit parameters (calculated from a least-squares 
fit of log 1 o(A r Jo' ( ! ) ) against logj^M^gQ*)) for the mass-richness re- 
lation |4j of simulated clusters with N^qq ^ 9 is compared to 
the best-fitting parameters for the SDSS maxBCG clusters with 
masses measured by |Johnston et al.| | |2007b} . 
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Table 5. The best-fit parameters (calculated from a least-squares 
fit of log 10 (AT|qq) against log 10 {M£oQ Ln )) for the mass-richness 
relation Q of simulated clusters with N^qq 9 is compared 
to the best-fit parameters for the SDSS maxBCG clusters with 
masses measured by |Reyes et al.| ( |2008[ > . 
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maxBCG cluster sample has been performed by |Reyes et al.| 
(2008). In Fig. [7| we compare the average cluster masses 
(^2uo an ) °f our simulated clusters to their results as a func- 
tion of richness. As for (M^oo), t ne average cluster masses 
(M^o an ) in the MS model are up to 30% larger than the 
SDSS maxBCG cluster masses, whereas the cluster masses 
of the WMAP3 models are comparable to those in the SDSS. 

The mass-richness relations shown in Fig. [6] and [7] sug- 
gest a power law, although with a steeper slope for (iVf 00 ) < 
10 than for (N%) > 10. Here, we fit the mean mass-richness 
relation for clusters with N^qq ^ 9 by: 



ti «-crit/mean / i» rgal \ * y-crit/mcan 

M 200 l JV 200/ — M 200I20 



N: 



200 

20 



crit/n 



(4) 



by Mandelbaum et al. ( 2008 I and Rozo et al. (2009a 



The best-fit parameters for our simulated catalogues are 
compared to those for the SDSS in Tables [4] and [5] 

The lower cluster masses in the WMAP3 models than 
in the MS and WMAP1-A models again reflect the different 
cosmologies. The less evolved dark matter structure in the 
WMAP3 cosmology requires more efficient star formation 
in order to match observed galaxy numbers. This results in 
more ridgeline galaxies in a dark matter halo of given mass. 
Consequently, for a given richness, halos are less massive in 
the WMAP 3 models than in the MS and WMAP1-A models. 

Except for the largest richness bin, observed cluster 
masses are smaller than in the MS model. This suggests 
a normalisation erg < 0.9 for our Universe (again consistent 
with the recent estimate as = 0.81 ± 0.03 by Komatsu et al. 
20091. 
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luminosity (£§ 00 / 



r gal v 



!,} of ridgeline galaxies within B|qq. Results for 
our MS and WMAP simulations are compared to SDSS results 
based on cluster masses by |Johnston et al.| p007b[ l. 



Table 6. The best-fit parameters for the mass-luminosity re- 
lation ^ in our simulations (calculated from a least-squares fit 
of Iogio(-^§nn) vs. log^A/jog )) are compared to the values by 
| Johnston et al.| ( |2007b| l for the SDSS maxBCG cluster sample. 
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3.4 Mass-luminosity relation 

The mass-luminosity relation computed by binning our sim- 
ulated cluster samples in luminosity L\ g is shown in Fig. [i] 
For given luminosity, the MS model yields cluster masses 
very similar to the SDSS masses of Johnston et al. (2007b]), 
whereas the WMAP3 model produces mean cluster masses 
that are generally smaller than the SDSS masses. The best- 
fit parameters for the mass-luminosity relation 



crit 
200 1 40 



real 

-^200 



4 x 10 n /i- 2 L 



(5) 



are listed in Table [6] 

The differences in the mass-luminosity relation between 
the various galaxy formation models can be explained in the 
same way as the differences in the mass-richness relation. 
The higher star formation efficiency in the WMAP3 models 
creates more bright ridgeline galaxies in a halo of a given 
mass A/So than in the MS and WMAP1-A models. This 
leads to lower average cluster masses at given luminosity 



3.5 The scatter in the mass-richness relation 

To compute a cluster mass function from cluster counts and 
mean masses as a function of richness, one needs to model 
the scatter in mass at each richness. In Fig. [9] distributions 
of the logarithm of cluster mass log 10 (M2oo) are shown for 
the MS model for several bins in N^qq- These distributions 
are well described by gaussians. 

For our various galaxy formation models, the standard 



deviation a 



1°E10( M 20()) 



of the scatter in the logarithm of clus- 



ter mass log 10 (M|oo ) a t gi ven N ^nn is listed for various y 



rgal 



and Lfpg-bins is in Tables 







14 



The scatter decreases with 



increasing iVfo or Lfoo ant ^ tends to be larger at given N^qq 
than at the corresponding Lfao- 

Our values for the scatter in the mass-richness rela- 
tion are in good agreement with those found empirically by 



Rozo et al. (2009a I for the maxBCG cluster sample (using 
X-ray luminosities as an additional mass proxy): We find 



'°8l0 «) 



0.15 



■ 0.2 for the richness bins with iVf^ 



while 



Rozo et al. 
i 



> 9, 
0.20 ± 0.09 



|2009a| find <7 logloW) 

for N$oo « 40. Moreover, our values are consistent with the 
scatter in the velocity dispersion-richness relation derived 



by |Becker et al. (20071 for the maxBCG clusters, if centre 
misidentification is taken into account. Note that there may 
be additional effects that increase the observed scatter, but 
are not modelled here. 



3.6 The mass function 

When the cluster abundance and the mass distribution at 
each richness is known, it is straightforward to reconstruct 
the cluster mass function ( Rozo et al.|2 009a). The differen- 
tial cluster number density (or differential mass function) is 
then given by a sum over all richness bins: 



dn(M 2 c £ t )/dM 2 c 5 t 



Y, mpdf, (MISS), (6) 



where rn denotes the space density and pdf ^(Mloo) the mass 
distribution for clusters in richness bin i. 

The results of Sec. 13.51 show that the mass distribu- 
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lustrates for the MS model, the reconstructed mass : 
matches the true mass function well for Mfoo > 2 x 10 14 . For 
smaller masses, the richness-selected cluster sample becomes 
incomplete in mass, and thus the reconstruction fails. 

The reconstructed cluster mass functions for the differ- 
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Figure 12. Comparison of the average weak-lensing mass profiles 
AS(i?) as a function of radius R for all clusters (solid line), for 
the correctly centred clusters (dashed line), and for the incorrectly 
centred clusters (dashed line) in the MS. 



ent galaxy models and the SDSS are compared in Fig. |11| 
The figure illustrates clearly why we expect a ACDM cos- 
mology with normalisation 0.72 < as < 0.9 to provide a 
better fit to the SDSS cluster data than the models consid- 
ered here. The values for the WMAP3 models are always 
much smaller those reconstructed from the SDSS, while the 
MS yields values above the observations. The reconstructed 
cluster mass function for the WMAP1 model generally fol- 
lows the MS results, but is visibly affected by sampling noise 
for larger cluster masses. 

Since the cluster mass function can be recovered from 
the cluster abundances and cluster masses as functions of 
richness, these quantities cannot vary independently if the 
cluster mass function is fixed. Different assumptions about 
the galaxy formation physics or the richness measurements 
that, for given richness, lead to higher cluster abundances 
will also yield lower cluster masses (and vice versa). Thus, 
the abundance-richness relation discussed in Sec. l3.2l and the 
mass-richness relation discussed in Sec. |3.3| provide comple- 
mentary information on the cosmology. 

3.7 Fits to the cluster mass profiles 

The results of Sec. |3.5| justify the use of a log-normal mass 



distribution for fits to observed mass profiles (e.g., by John- 
ston et al.|2007b||Reyes et al.|2008| ). Here, we illustrate tha ; 
one can indeed obtain a good fit to the simulated mean mass 
profiles of clusters with assumptions similar to those used, 
e.g., by Johnston et al. (2007b I. 

The mass profiles discussed in Sec. |3.1| are a mixture 
of correctly and incorrectly centred clusters. In Fig. |12| the 
simulated average mass profiles of all clusters are compared 
to the profiles of correctly and incorrectly centred clusters. 
The profiles agree well for large radii, but differ significantly 
below a certain radius R w 0.5/i -1 Mpc for clusters with 
= 3, and R w l/i -1 Mpc for clusters with 



richness A|qq 



As Fig. 1 12| illustrates, centre misidentification has a sig- 
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Figure 13. Fit to the weak-lensing mass profile AS(iJ) as a 
function of radius R of correctly centred clusters in the MS model. 
Shown are the measured profiles (solid line), the 3-component fit 
(dashed line), the central galaxy contribution (dash dotted line), 
the DM halo contribution (dotted line) , and the contribution from 
neighbouring masses (dash-dot-dotted line). 



nificant impact on the average cluster mass profiles and thus 
needs to be taken into account in profile fits. Since stack- 
ing weak-lensing mass profiles is linear, we can discuss the 
contributions from correctly and incorrectly centred clusters 
separately. A weighted average of fits to these two compo- 
nents constitutes a fit to the average mass profile of all clus- 
ters in a richness bin. 

We assume that the average mass profile of the correctly 
centred clusters in a richness bin consists of a central galaxy 
component, which we model as point mass, a mean dark 
matter halo modelled as an average over spherical NFW 



profiles ( Navarro et al. 1997 1, and a contribution from neigh- 
bouring mass concentrations. A log-normal distribution with 
mean (Affoo) an d standard deviation o"i ogi0 (M crit ) gi ven by 
Table [7] is assumed for the halo masses. Furthermore, we 
assume that the concentration c of halos with mass Mfoo 
follows a log-normal distribution with mean 



(c) (Mwo) 



4.67 



lO^/i^M© 



(7) 



and standard deviation o"i ogl0 ( c ) = 0.15 ( Neto et al.||2007 l 



The resulting fit to the average weak-lensing mass pro- 
file of correctly centred clusters is shown in Fig. [13] for clus- 
ters in the MS model with richness 18 sC Afoo < 25. A fit 
of similar quality can also be obtained for the other richness 
bins. 

We now turn to the contribution from incorrectly cen- 
tred clusters. Figure [12] shows that the average profile of the 
incorrectly centred clusters increases with decreasing radius 
even for small radii. This is a consequence of the particu- 
lar choice for the apparent centre of these clusters in our 
simulation, namely another massive cluster galaxy within a 
massive subhalo. If the apparent cluster centres were cho- 
sen randomly, one would expect the profile to decrease with 
decreasing radius for R < l/i -1 Mpc ( Johnston et al. 2007b I. 




0.1 1 

R [h~ l Mpc] 

Figure 14. Fit to the weak-lensing mass profile AE(_R) as a func- 
tion of radius R of incorrectly centred clusters in the MS model. 
Shown are the measured profiles (solid line), the 4-component fit 
(short-dashed line), the central galaxy contribution (dash-dotted 
line), the DM halo contribution (dotted line), the contribution 
from neighbouring masses (dash-dot-dottcd line), and the sub- 
halo contribution (long-dashed line). 



correctly centred clusters, we thus need four components: a 
central galaxy component (again modelled as point mass), 
a mean dark matter halo modelled as an NFW profile con- 
volved with a 2D Gaussian, a contribution from neighbour- 
ing masses, and a subhalo, which we model as a truncated 
NFW profile ( jBaltz et alJ]2009| ). As Fig. 



14 



illustrates, the 

subhalo component is essential for a good fit to the simu- 
lated mass profiles. 



4 SUMMARY AND DISCUSSION 

In this work, we have used A r -body simulations of cos- 
mic structure formation with semi-analytic galaxy forma- 
tion modelling to test this modelling and to investigate how 
the properties of optically selected galaxy groups and clus- 
ters depend on cosmological parameters. We have created 
catalogues of simulated galaxy groups/clusters from model 
galaxy catalogues (by De Lucia and Blaizot 2007 Wang 
et al.||2008 1. We have computed weak-lensing mass profiles 
and various other properties for these clusters as a function 



of cluster richness iVfgg and luminosity L|qq (as defined by 



gal 



the maxBCG algorithm of Koester et al.||2007b[ ), and com- 
pared the results to observations of clusters in the SDSS 
||Sheldon et al.||2009"a| | Johnston et al.|[2007b"] |Reyes et ah 
20081) 



To obtain a good fit to the average mass profiles of in- 



We find that the simulated weak-lensing mass profiles 
and the observed profiles of the SDSS maxBCG clusters 
agree remarkably well in detailed shape and amplitude. 
Moreover, all simulations reproduce the observed mass- 
richness relation within ~ 30% over the whole range probed 
by the SDSS clusters. The MS and the WMAP1-A sim- 
ulation also yield cluster abundances very similar to the 
observed abundances. This shows that the models consid- 
ered here provide a good description of the masses, den- 
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sity profiles and optical properties of galaxy clusters as well 
as the optical cluster selection and richness estimation of 
the maxBCG algorithm. Evidently, using mock galaxy cata- 
logues based on a large high-resolution ACDM structure for- 
mation simulation and semi-analytic galaxy formation mod- 
els makes it possible to create very realistic mock cluster 
catalogues for surveys like the SDSS (when problems with 
the ridgeline colours are overcome). 

Although the underlying A r -body simulations assume 
different cosmological parameters that lead to different DM 



halo abundances (Springel et al. 2005 Wang et al. 20081, 



all the galaxy models used here are able to reproduce the 
observed abundance and two-point correlations of galaxies 
reasonably well through careful adjustment of their star- 



formation efficiency and feedback parameters ( De Lucia and 



Blaizot||2007| |Wang et al.||2008[ ). They differ, however, in 
their predicted cluster abundance as a function of iVfoo- 
The MS model and the WMAP1-A model, both of which 
use cosmological parameters based on lst-year WMAP data 



(Spergel et al. 20031, predict cluster abundances that are 
compatible with the observed values ( Sheldon et al.|2009a l, 
whereas the WMAP3-A and WMAP3-B models, which use 
parameters consistent with the WMAP 3rd-year results 



(Spergel et al. 20071, yield abundances that are lower by 



a factor 2-3. 

The cluster masses predicted as a function of richness 
Afoo or luminosity L^qq also differ for the different galaxy 
formation models. At given richness or luminosity, the MS 
and WMAP 1- A models produce clusters that are up to 
30% more massive than observed, while the WMAP3 cluster 
masses are similar to or lower than the observed masses. 

The different abundances and average cluster masses 
in our various galaxy formation models are primarily a re- 
flection of the different underlying cosmologies. Because ha- 
los in the WMAP3 cosmology are less massive than in the 
MS/WMAP1 cosmology, the WMAP 3 models need more ef- 
ficient star formation to match the observed galaxy number 
densities. This produces more ridgeline galaxies in a dark 
matter halo of given mass. Consequently, clusters at a given 
richness are less massive in the WMAP3 models than in the 
MS and WMAP1-A models. Nevertheless, the higher star 
formation efficiency in the WMAP3 models does not fully 
compensate for the lower number of massive halos. So there 
are fewer rich clusters in the WMAP3 models than in the 
MS and WMAP 1- A models. 

The lower cluster masses in the observations compared 
to the MS model suggests that our Universe would be bet- 
ter described by as < 0.9. The higher observed rich clus- 
ter abundance than in the WMAP3 models suggests that 
as > 0.72. Thus, more precise agreement between predicted 
and observed cluster properties is expected for an interme- 
diate value 0.72 < as < 0.9. This corroborates the findings 
by |Rozo et"aT] ( |2009b[ ) that the SDSS maxBCG cluster data 
favour 08 ~ 0.83. 

Our results confirm that the mass distribution of clus- 
ters of given richness is well described by a log-normal dis- 
tribution. This justifies both the assumption of such distri- 
butions and the specific scatter values adopted in previous 
work which modelled stacked cluster mass profiles or recon- 
structed cluster mass functions. Fits to the stacked mass pro- 
file of clusters whose centre has erroneously been identified 
with a non-central cluster galaxy should take into account a 



halo component associated with this non-central galaxy in 
addition to its stellar mass, the halo component of the main 
cluster, and surrounding large-scale structure. 

Our simulations required many simplifying assumptions 
about the richness measurements (e.g. about limiting mag- 
nitudes, colours, or projection effects), which could result 
in biased estimates. Although predictions for cluster abun- 
dances or for cluster masses are individually subject to such 
modelling errors, they do not vary independently for varying 
assumptions about the richness measurement. For example, 
assuming a fainter magnitude limit results in lower cluster 
masses and higher cluster numbers. Thus, it proves difficult 
to decrease average cluster masses for MS model in order to 
better match the observations without producing too many 
rich clusters. Similarly, the cluster numbers in the WMAP3 
models can only be brought into agreement with the ob- 
served counts at the cost of cluster masses which are too 
small. 

Similar reasoning reveals that changes to the galaxy 
formation description alone, though they could change the 
number, colour, or brightness of the model galaxies, could 
only lead to better agreement for either cluster abundance or 
cluster mass, but not both. This stems from the constraint 
that the abundances and masses of clusters as a function 
of richness must 'add up' to reproduce the underlying clus- 
ter mass function regardless of the specific galaxy formation 
model. To reach better agreement for both, the number of 
massive clusters has to be adjusted, too, by changing the 
cosmology. 

Our results demonstrate that, on the one hand, the 
cluster mass-richness relation and the cluster abundance- 
richness function together provide strong constraints for the 
cosmology even without perfect knowledge of the galaxy for- 
mation physics. On the other hand, our findings show that 
the cluster abundance and cluster mass as functions of rich- 
ness can also be used directly, in addition to galaxy abun- 
dance and galaxy two-point correlation, to test the galaxy 
models. 

In future work, one should test how much the agree- 
ment between observed and predicted cluster properties can 
be improved by choosing different cosmological parameters 
for the simulations (e.g. the values currently favoured by 
other observations such as Komatsu et al. 2009). In addi- 



tion, the galaxy models discussed here should be improved 
to better match the observed colours of ridgeline galaxies at 
higher redshift. (If we had not adjusted the ridgeline colour 
selection by hand, the models would have contained almost 
no clusters with z > 0.25.) Future simulations should also 
probe other cosmologies and galaxy models. These simula- 
tions need to have a much larger volume than the WMAP 
simulations to have good enough statistics to match up- 
coming observations. (The statistical errors on the mass- 
richness relation in the WMAP simulations are comparable 
to the uncertainties in current observations.) 

More realistic modelling of cluster selection and char- 
acterisation could be achieved by running the observational 
cluster-finding algorithms on mock galaxy catalogues cre- 
ated from the simulations. Moreover, ray-tracing techniques 
could be used to simulate realistically the weak-lensing mass 
measurements and to assess their statistical accuracy and 
possible systematic uncertainties. 

Finally we note that there is some tension between our 
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finding that the abundances and weak-lensing masses of 
clusters favour cosmologies with a normalisation as ~ 0.8 
over those with og = 0.72 (in agreement with estimates by 



Rozo et al. 



by Li et al. 



2009b Komatsu et al. 20091 and the findings 



(2009) and Cacciato et al. (20091 that galaxy- 



galaxy lensing and galaxy clustering data are consistent with 
as = 0.73. It is beyond the scope of this paper to analyse 
possible reasons for this discrepancy and whether it can be 
resolved with better handling of measurement systematics 
or improved structure formation models, but this should be 
done in future work. 
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Table 7. The average properties of the model clusters in the MS as inferred from the galaxy model by |De Lucia and B laizot ( 2007J) in 
bins of richness N^qq: the comoving number density n, the average number {Nf^) of ridgeline galaxies, the average number (N^qq) of 

ridgeline galaxies within fifoo» ^ ne avera g e cluster radius (-Rfoo)' *' le average ridgeline i-band luminosity (Lfoo) within fifoO' * ne avera g e 
cluster masses {M%qq), and the standard deviation cr logl j M critj of the logarithmic cluster mass. For bins containing < 7 independent 
clusters, <r, /.raits is omitted. 



bin JVgJ n [^Mpc-3] (N.^) (AgJ) (fl**) [fc-iMpc] [ft-lL Q ] <Mf$> [ft-iM ] ^ logl0 ( M «jt) 
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Table 8. The properties of the cluster sample in the WMAP1-A model in bins of richness JVfgg (see Table [jj for a description of the 
listed quantities). 
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Table 9. The properties of the cluster sample in the WMAP3-B model in bins of richness N^qq (see Table [jj for a description of the 
listed quantities). 
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Table 10. The properties of the cluster sample in the WMAP3-C model in bins of richness N§qq ( see Table [7] for a description of the 
listed quantities). 
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8 


5.6 x IO" 6 


7.3 


8.0 


0.68 


1.7 x 10 11 


2.9 x 10 13 


0.24 


7 


9-11 


7.9 x IO" 6 


9.3 


9.7 


0.77 


2.1 x 10 11 


3.9 x 10 13 


0.24 


8 


12-17 


6.6 x IO" 6 


14.6 


14.2 


0.97 


3.1 x 10 11 


6.8 x 10 13 


0.19 


9 


18-25 


2.6 x IO" 6 


20.4 


20.5 


1.16 


4.3 x 10 11 


9.5 x 10 13 


0.20 


10 


26-40 


1.2 x 10~ 6 


30.1 


29.6 


1.41 


6.5 x 10 11 


1.5 x 10 14 




11 


41-70 


2.5 x 10~ 7 


65.7 


53.2 


1.77 


1.1 x 10 12 


2.5 x 10 14 




12 


71-220 


2.7 x IO" 7 


77.5 


82.0 


2.09 


1.6 x 10 12 


3.7 x 10 14 





Table 11. The average properties of the model clusters in the MS model in bins of luminosity £|qq (see Table [jj for a description of 
the listed quantities). 



bin Lf a > [10 1 °h- 2 L Q ] n [^Mpc" 3 ] (N^) (iV|g*> <Rf^> [h^Mpc] (Lf&) [h-*L ] (M^) [HMg] tr } 



1 


5.00 - 6.24 


4.1 


X 


10~ 5 


2.9 


3.3 


0.39 


5.6 


X 


10 10 


1.0 


X 


10 13 


0.27 


2 


5.24 - 7.80 


4.4 


X 


10~ 5 


3.3 


3.6 


0.42 


7.0 


X 


10 10 


1.3 


X 


10 13 


0.27 


3 


7.80 - 9.74 


3.9 


x 


10~ 5 


3.8 


4.1 


0.46 


8.7 


X 


10 10 


1.8 


X 


10 13 


0.27 


1 


9.74 - 12.2 


3.1 


x 


10~ 5 


4.6 


4.9 


0.51 


1.1 


X 


10 11 


2.4 


X 


10 13 


0.26 


5 


12.2 - 15.2 


2.2 


x 


10~ 5 


5.7 


5.9 


0.58 


1.4 


X 


10 11 


3.2 


X 


10 13 


0.24 


6 


15.2 - 19.0 


1.6 


X 


10~ 5 


7.2 


7.1 


0.66 


1.7 


X 


10 11 


4.2 


X 


10 13 


0.21 


7 


19.0 - 23.7 


1.1 


X 


10" 5 


9.2 


9.3 


0.76 


2.1 


X 


10 11 


5.6 


X 


10 13 


0.19 


8 


27.7 - 29.6 


7.7 


X 


10~ 6 


11.7 


11.8 


0.87 


2.6 


X 


10 11 


7.1 


X 


10 13 


0.18 


9 


29.6 - 36.9 


5.2 


X 


10~ 6 


14.9 


15.0 


0.99 


3.3 


X 


10 11 


9.5 


X 


10 13 


0.16 


10 


36.9 - 46.1 


3.7 


X 


10" 6 


18.7 


18.9 


1.12 


4.1 


X 


10 11 


1.2 


X 


10 14 


0.16 


11 


46.1 - 57.6 


2.5 


X 


10~ 6 


23.7 


23.8 


1.26 


5.1 


X 


10 11 


1.6 


X 


10 14 


0.15 


12 


57.6 - 71.9 


1.6 


X 


10~ 6 


29.3 


29.8 


1.40 


6.4 


X 


10 11 


2.0 


X 


10 14 


0.13 


13 


71.9 - 89.8 


9.6 


X 


io- 7 


37.0 


37.9 


1.57 


8.0 


X 


10 11 


2.6 


X 


10 14 


0.13 


11 


89.8 - 112.1 


6.5 


X 


io- 7 


46.7 


47.5 


1.76 


1.0 


X 


10 12 


3.3 


X 


10 14 


0.12 


15 


112.1 - 140 


4.1 


X 


IO" 7 


56.8 


59.7 


1.94 


1.2 


X 


10 12 


4.1 


X 


10 14 


0.13 


16 


140 - 450 


3.6 


X 


io- 7 


84.7 


91.2 


2.31 


1.9 


X 


10 12 


6.4 


X 


10 14 


0.16 
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Table 12. The average properties of the model clusters in the WMAP1-A model in bins of luminosity i|go ( see Table|jJfor a description 
of the listed quantities). 



bin Lf a > [l O l°fc- 2 L ] n [^Mpc" 3 ] (N*£) (iVf^) (Rf^> [fc-iMpc] <Lf^> [h^L e ] (M-*) [k^Mg] a, 



1 


5.00 - 6.24 


4.4 


X 


10~ 5 


2.9 


3.3 


0.39 


5.6 


X 


10 10 


9.8 


X 


10 12 


0.27 


2 


5.24 - 7.80 


5.0 


X 


10~ 5 


3.2 


3.6 


0.42 


7.0 


X 


10 10 


1.3 


X 


10 13 


0.26 


3 


7.80 - 9.74 


4.4 


X 


10~ 5 


3.9 


1.2 


0.46 


8.7 


X 


10 10 


1.7 


X 


10 13 


0.26 


1 


9.74 - 12.2 


3.7 


X 


10~ 5 


4.7 


5.0 


0.52 


1.1 


X 


10 11 


2.2 


X 


10 13 


0.24 


5 


12.2 - 15.2 


2.5 


X 


10~ 5 


5.6 


5.9 


0.57 


1.4 


X 


10 11 


2.8 


X 


10 13 


0.22 


6 


15.2 - 19.0 


1.7 


X 


10~ 5 


7.2 


7.2 


0.65 


1.7 


X 


10 11 


3.8 


X 


10 13 


0.21 


7 


19.0 - 23.7 


1.2 


X 


10~ 5 


9.1 


9.2 


0.75 


2.1 


X 


10 11 


5.0 


X 


10 13 


0.19 


8 


27.7 - 29.6 


9.3 


X 


IO" 6 


11.4 


11.5 


0.85 


2.6 


X 


10 11 


6.6 


X 


10 13 


0.17 


9 


29.6 - 36.9 


5.9 


X 


10~ 6 


14.9 


14.5 


0.97 


3.3 


X 


10 11 


8.5 


X 


10 13 


0.17 


10 


36.9 - 46.1 


4.3 


X 


10~ 6 


21.0 


19.8 


1.15 


4.1 


X 


10 11 


1.1 


X 


10 14 


0.15 


11 


46.1 - 57.6 


2.5 


X 


10~ 6 


22.8 


24.0 


1.26 


5.1 


X 


10 11 


1.3 


X 


10 14 




12 


57.6 - 71.9 


1.4 


X 


IO" 6 


29.2 


30.4 


1.42 


6.5 


X 


10 11 


1.8 


X 


10 14 




13 


71.9 - 89.8 


1.7 


X 


10~ 6 


36.2 


37.6 


1.59 


8.0 


X 


10 11 


2.2 


X 


IO 14 




14 


89.8 - 112.1 


3.9 


X 


io- 7 


39.4 


43.5 


1.65 


9.6 


X 


10 11 


2.3 


X 


IO 14 




15 


112.1 - 140 


1.8 


X 


10- 7 


55.0 


63.9 


1.88 


1.2 


X 


10 12 


2.3 


X 


IO 14 




16 


140 - 450 


5.1 


X 


io- 7 


94.5 


108.1 


2.50 


2.1 


X 


10 12 


1.0 


X 


10 15 





Table 13. The average properties of the model clusters in the WMAP3-B model in bins of luminosity L| q (see Tablejjjfor a description 
of the listed quantities). 



bin LS> [10 10 / l - 2 L Q ] n [^Mpc" 3 ] (N*£) (jVf a J> (H**) [h^Mpc] (L^ Q ) [h^L e ] (M%%) [h^U @ ] a, 



1 


5.00 - 6.24 


3 


.0 


X 


10~ 5 


2.8 


3.3 





.39 


5.6 


X 


10 10 


8 


.5 


X 


10 12 


0.29 


2 


5.24 - 7.80 


3 


.2 


X 


10~ 5 


3.2 


3.6 





.41 


7.0 


X 


10 10 


1 


.1 


X 


10 13 


0.29 


3 


7.80 - 9.74 


3 


.0 


X 


10~ 5 


3.8 


4.0 





.45 


8.7 


X 


10 10 


1 


.5 


X 


10 13 


0.29 


1 


9.74 - 12.2 


2 


.1 


X 


IO" 5 


1.4 


1.6 





.49 


1.1 


X 


10 11 


1 


.9 


X 


10 13 


0.27 


5 


12.2 - 15.2 


1 


.4 


X 


IO" 5 


5.3 


5.4 





.54 


1.4 


X 


10 11 


2 


.4 


X 


10 13 


0.26 


6 


15.2 - 19.0 


8 


.9 


X 


10~ 6 


6.8 


6.8 





.63 


1.7 


X 


10 11 


3 


.3 


X 


10 13 


0.24 


7 


19.0 - 23.7 


5 


.1 


X 


10~ 6 


8.9 


8.7 





.74 


2.1 


X 


10 11 


4 


.2 


X 


10 13 


0.22 


8 


27.7 - 29.6 


3 


.7 


X 


IO" 6 


11.2 


10.8 





.83 


2.6 


X 


10 11 


5 


.9 


X 


10 13 


0.22 


9 


29.6 - 36.9 


2 


.6 


X 


IO" 6 


14.0 


13.6 





.93 


3.3 


X 


10 11 


7 


.5 


X 


10 13 




10 


36.9 - 46.1 


1 


.8 


X 


10"° 


18.4 


17.5 


1 


.06 


4.1 


X 


10 11 


1 


.0 


X 


10 14 




11 


46.1 - 57.6 


9 


.6 


X 


io- 7 


25.5 


22.2 


1 


.22 


5.1 


X 


10 11 


1 


.4 


X 


10 14 




12 


57.6 - 71.9 


5 


.0 


X 


io- 7 


30.2 


26.2 


1 


.32 


6.4 


X 


10 11 


1 


.6 


X 


10 14 




13 


71.9 - 89.8 


2 


.6 


X 


io- 7 


43.8 


35.4 


1 


.15 


7.9 


X 


10 11 


2 


.1 


X 


IO 14 




11 


89.8 - 112.1 


1 


.5 


X 


io- 7 


59.7 


49.1 


1 


.68 


1.0 


X 


10 12 


2 


.3 


X 


IO 14 




15 


112.1 - 140 


1 


.3 


X 


io- 7 


67.9 


61.4 


1 


.85 


1.2 


X 


10 12 


2 


.9 


X 


IO 14 




16 


140 - 450 


1 


.4 


X 


HT 7 


72.6 


81.3 


2 


.03 


1.5 


X 


10 12 


3 


.8 


X 


IO 14 
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Table 14. The average properties of the model clusters in the WMAP3-C model in bins of luminosity £§qq (see Tablej^jfor a description 
of the listed quantities). 



bin Lf»> [l O l°fc- 2 L ] n [^Mpc" 3 ] (N^) (iV|gJ> [fc^Mpc] (L*g,> [ft-iL ] (M-") [fc- 1 ^] ^o gl0 (M| 5 «) 



1 


5.00 - 6.24 


3.9 


X 


10~ 5 


2.9 


3.3 


0.39 


5.6 


X 


10 10 


8.1 


X 


10 12 


0.29 


2 


5.24 - 7.80 


3.9 


X 


10~ 5 


3.3 


3.6 


0.42 


7.0 


X 


10 10 


1.0 


X 


10 13 


0.30 


3 


7.80 - 9.74 


3.3 


X 


10~ 5 


3.8 


1.1 


0.46 


8.7 


X 


10 10 


1.3 


X 


10 13 


0.29 


1 


9.74 - 12.2 


2.6 


X 


10~ 5 


4.5 


4.8 


0.50 


1.1 


X 


10 11 


1.7 


X 


10 13 


0.27 


5 


12.2 - 15.2 


1.6 


X 


IO" 5 


5.5 


5.9 


0.56 


1.4 


X 


10 11 


2.2 


X 


10 13 


0.24 


6 


15.2 - 19.0 


1.0 


X 


10~ 5 


6.7 


7.0 


0.63 


1.7 


X 


10 11 


2.9 


X 


10 13 


0.23 


7 


19.0 - 23.7 


6.2 


X 


10~ 6 


8.8 


8.8 


0.73 


2.1 


X 


10 11 


3.8 


X 


10 13 


0.21 


8 


27.7 - 29.6 


4.7 


X 


IO" 6 


11.6 


11.7 


0.86 


2.6 


X 


10 11 


5.3 


X 


10 13 


0.20 


9 


29.6 - 36.9 


3.5 


X 


10~ 6 


14.9 


15.0 


0.99 


3.3 


X 


10 11 


6.9 


X 


10 13 


0.17 


10 


36.9 - 46.1 


2.1 


X 


10~ 6 


18.3 


17.9 


1.09 


1.1 


X 


10 11 


9.3 


X 


10 13 




11 


46.1 - 57.6 


1.1 


X 


10~ 6 


23.3 


23.2 


1.24 


5.1 


X 


10 11 


1.3 


X 


10 14 




12 


57.6 - 71.9 


6.4 


X 


io- 7 


31.2 


27.9 


1.38 


6.4 


X 


10 11 


1.5 


X 


10 14 




13 


71.9 - 89.8 


3.7 


X 


io~ 7 


37.9 


33.9 


1.49 


7.8 


X 


10 11 


1.8 


X 


10 14 




14 


89.8 - 112.1 


1.4 


X 


io- 7 


57.8 


48.0 


1.70 


1.0 


X 


10 12 


2.2 


X 


10 14 




15 


112.1 - 140 


1.1 


X 


io- 7 


67.5 


59.7 


1.87 


1.2 


X 


10 12 


2.7 


X 


10 14 




16 


140 - 450 


2.8 


X 


io- 7 


77.3 


81.1 


2.07 


1.6 


X 


10 12 


3.7 


X 


IO 14 





